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We present an explicit analytic calculation of the energy-level shift of an atom in front of a non- 
dispersive and non-dissipative dielectric slab. We work with the fully quantized electromagnetic 
field, taking retardation into account. We give the shift as a two-dimensional integral and use 
asymptotic analysis to find expressions for it in various retarded and non-retarded limiting cases. 
The results can be used to estimate the energy shift of an atom close to layered microstructures. 
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I. INTRODUCTION 

Control and manipulation of cold atoms have become 
fundamentally important due to their central role in the 
development of nanotechnology and as a tool for inves- 
tigating the mechanisms underlying macroscopic mani- 
festations of quantum physics It now seems feasi- 
ble to control them on a /im length-scale by utilizing 
microstructured surfaces — also known as atoms chips 
— with promising areas of application such as quantum 
information processing with neutral atoms, integrated 
atom optics, precision force sensing, and studies of the 
interaction between atoms and surfaces. 

For this reason, e.g., experiments using Bose-Einstein 
condensates for measuring the Casimir-Polder force [2( 
have been developed. Typically, the dielectric substrate 
utilized in such experiments [3| carries a very thin top 
layer of another material, generally graphite or gold. 
However, in explicit analytic theories, this finite thick- 
ness has often been neglected and the system has been 
treated as a semi-infinite half-space. Here we are aiming 
at an approach that lets us include the thickness of such 
a layer as a parameter into the calculation, allowing us to 
obtain analytic expressions for the Casimir-Polder force 
on an atom. 

We are going to consider a ground-state atom close to 
a non-dispersive dielectric slab, which is one of the fews 
systems of high symmetry for that the electromagnetic 
field can be quantized through an exact normal-mode 
expansion with manageable effort 0]. The assumption 
of absent dispersion and absorption is working well for 
all but very few systems, namely those where the atom 
has a strong transition very near an absorption line in 
the medium, which in practice is something very diffi- 
cult to engineer. Nevertheless, if one wishes to consider 
more complex systems with atoms near absorbing bound- 
aries, the calculation will require other methods to study 
quantum electrodynamics [5(. However, applying such 
results to a particular case might require extensive nu- 
merical calculations for the case at hand. Using those 
types of formalism, the energy-level shift of atoms due to 
the presence of media with diverse magneto-electric prop- 
erties has been calculated for several systems Q. Here, 
by contrast, we are not aiming at general expressions. 



The focal point of the present work is to obtain simple 
and practical formulae that are useful for estimates and 
can be applied very easily to experimental situations. 

The energy shift in an atom close to a dielectric slab 
comes about due to its interaction with electromagnetic 
field fluctuations, which in turn are affected by the pres- 
ence of the slab. Thus, a quantization of the electro- 
magnetic field in the presence of a layered system is re- 
quired. Even though the field quantization for this sys- 
tem had been studied previously we have recently 
re-considered the problem and provided a proof of the 
completeness of the electromagnetic field modes that was 
missing in previous works on the problem. This proof of 
completeness is very useful in that it removes any am- 
biguity in how to normalize and sum over the electro- 
magnetic field modes, and in this way also establishes 
the correct density of states which had previously been a 
subject of disagreement [8j|. 

By solving the Helmholtz equation and imposing the 
corresponding continuity conditions at the faces of the 
slab, it was shown in [4], 0| that the field modes for this 
system comprise of travelling and trapped modes. The 
travelling modes have a continuous frequency spectrum, 
and are composed of incident, reflected and transmitted 
parts outside the slab. The trapped modes arise due to 
solutions of the Helmholtz equation with purely imag- 
inary normal wave vector outside the slab. Physically, 
they come about due to repeated total internal reflection 
inside the dielectric, and emerge as evanescent fields out- 
side the slab. They exist only at certain discrete frequen- 
cies which depend on polarization direction and parity 
and are obtained through the dispersion relations. 

The atomic energy shift is obtained by means of 
second-order perturbation theory, and involves a product 
of electromagnetic mode functions which is summed over 
intermediate virtual photon states. Thus the atomic en- 
ergy shift receives two quite separate contributions: one 
from the continuous set of travelling modes, and the other 
from the discrete set of trapped modes. The first is an in- 
tegral over wave numbers, and the second a sum over dis- 
crete wave numbers that satisfy a quite complicated dis- 
persion relation. In practice, both must be considered to- 
gether at all times to avoid divergent terms appearing in 
each separate contribution but cancelling between them. 
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FIG. 1: The geometry of the dielectric slab. 



This technically seemingly hopeless task can, however, be 
dealt with simply and elegantly by using the summation 
method of Q, which re-expresses both the integral over 
continuous modes and the sum over discrete modes as 
a single contour integral in the complex plane. Wc shall 
show below that this trick yields a closed-form expression 
for the atomic energy shift and permits easy asymptotic 
analysis for various regimes, yielding the kinds of simple 
formulae that we are after for estimating the effect of the 
layer thickness in experimental situations. 



II. DESCRIPTION OF THE SYSTEM 

We consider a dielectric slab of finite thickness L sur- 
rounded by vacuum, as is shown in Fig. [TJ We assume 
the material to be a non-dispersive and non-absorbing 
dielectric, which is a simple but good model for an im- 
perfectly reflecting material. Thus the material is charac- 
terized solely by its refractive index n, which is real and 
the same for all frequencies. While any real material has 
of course to be transparent at infinite frequencies, this 
non-dispersive model captures the essential properties of 
an imperfect reflector. In particular, it includes evanes- 
cent waves, whose absence in perfect-reflector models can 
be problematic @. 

Since the dielectric is homogeneous in the x and y di- 
rections, the dielectric permittivity of the configuration 
depends only on the z coordinate and is given by 



e(z) 



n 2 for - L/2 < z < L/2 

1 for \z\ > L/2 . 



We assume the atom to be neutral and in its ground 
state. Also, we shall make use of the electric dipole ap- 
proximation, which is adequate because, for the relevant 
modes, the electromagnetic field varies slowly over the 
size of the atom. We assume that the atom's center is 
fixed at the position r = (0, 0, z ). 

The model assumes that the interaction between the 
atom and the surface is purely electromagnetic, i.e. that 
there is negligible wave-function overlap between the 
atomic electron and the surface. 



We shall work with an interaction Hamiltonian be- 
tween the atom and the quantized electromagnetic field 
that is given by 



H int = -M'E(r,i) , 



(1) 



which is the lowest-order multipole Hamiltonian and cor- 
responds to the electric-dipole interaction. In this equa- 
tion, fi = e(r — ro) is the electric-dipole moment of 
the atomic electron, and E(r, t) is the transverse electric 
field. Unlike the minimal-coupling Hamiltonian p • A, 
the Hamiltonian |T|) includes the electrostatic interac- 
tion between the atomic dipole and its images on the 
other side of vacuum-dielectric interfaces. As shown pre- 
viously in a similar context fl"ot | , the Hamiltonian {T]) may 
be more convenient for calculations that aim to derive 
energy shifts in cases where the retardation of the elec- 
tromagnetic interaction matters. 

The quantization of the electromagnetic field has been 
discussed in detail previously [1, @], and thus, we shall 
only sketch the procedure here. We work with the elec- 
tromagnetic potentials $(r, t) and A(r, t) and choose the 
generalized Coulomb gauge 



V ■ [e(z)A(r)] = 0. 



(2) 



Furthermore, since the overall system is neutral, we can 
set <J>(r,i) = 0. Thus the field equations reduce to the 
wave equation for A(r, t) everywhere except right on 
the interfaces z — ±L/2. At the interfaces we solve 
Maxwell's equations directly by imposing the continuity 
conditions 



E h D x , B 



continuous . 



(3) 



In this way the electromagnetic field modes can be writ- 
ten, for travelling modes f„ R as left- or right-incident 
waves made up of incoming, reflected and transmitted 
parts, and for trapped modes fj?' A , as symmetric or an- 
tisymmetric waves inside the slab with evanescent fields 
outside. We list these modes in Appendix [Al 

Equipped with a complete set of solutions to the clas- 
sical field equations, we can proceed to quantize the elec- 
tromagnetic field by using the technique of canonical 
quantization, i.e. by introducing annihilation and cre- 
ation operators a„, a\,. Then the expansion for the elec- 
tric field operator E(r, t) in terms of the normal modes 
fj,(r) reads 

E(r,i) = i V f^L[a u e-^%{r) - bJ^CW] • 

(4) 

where the subscript v = (k, A) is a composite label in- 
cluding both the polarization A = TE, TM and the wave 
vector k. 



III. ENERGY-LEVEL SHIFT 

Since the interaction Hamiltonian ([T]) is linear in the 
electric field, whose vacuum expectation value vanishes, 
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there is no energy shift to first order in H[ nt . Therefore, 
the lowest-order contribution to the shift comes from 
second-order perturbation theory, so that the shift is of 
first order in the fine-structure constant a II ill. 



j;i/|/i-E(r,*)|*;0 



Ej — Ei 



(5) 



where E(r, t) is given by Eq. ((4}. In this equation, the in- 
termediate state \j; v) is a composite state with an atom 
in the excited state \j) and the electromagnetic field car- 
rying a photon of energy oj v . Similarly, the initial state 
\i; 0) describes an atom in its ground state \i) and the 
electromagnetic field in the vacuum state. In the electric 
dipole approximation we can write 



(j|f(r)-/x|i)~f(ro)-(iHi) 



(6) 



since the field varies slowly over the size of the atom and 
we can therefore assume that across the atom it is almost 
the same as at its center r = (0, 0, z ). With this, the 
shift reads 



AE = -—YY 



E, 



|C(ro)-OV|i>l . (7) 



where we have introduced the abbreviation Eji = Ej — 
Ei. Since \i) is a state of definite angular momentum, 
different components of fi lead to different intermediate 
states \j) that are mutually orthogonal, and the shift 
simplifies to 



AE = 



E 



|f*(r )| 2 |0>N) 



(8) 



Furthermore, we are going to abbreviate the moduli 
squares of the matrix elements of the dipole-momentum 
operator fi between the initial state i and the intermedi- 
ate states j, and distinguish only the components parallel 
and perpendicular to the slab, 



kri 2 = ioi/a>r 

lM||| 2 = IOK|i>l 2 



with a — x,y, z 

|2 



2 = I0KK>I 2 



The sum over v in Eq. ([5]) is a sum over all field modes, 
which, as explained earlier and easily seen from Appendix 
|A"1 comprise a continuous set of travelling modes and a 
discrete set of trapped modes. The contribution from the 
travelling modes gives rise to the shift 



A£ trav = -— V V 
2e ^ ^ 



£ 



jjii \=TE,TM <J 
R 



d J k- 



Eji 



x(|f^(ro)| 2 + |f£ ttr (r )| a ) |M, 



(9) 



where the sum over a runs over the x, y and z com- 
ponents of the dipole moment and of the polarization 



vector that is incorporated in the mode functions. As we 
are interested in the change in the energy levels of the 
atom solely due to the presence of the dielectric slab, we 
renormalize the energy shift and remove from Eq. §§§ 
the part that arises due to the interaction between the 
atom and the electromagnetic field in free space, i.e. the 
Lamb shift. Conveniently, this rids the calculation of 
any divergences, provided travelling and trapped modes 
are considered together (cf. e.g. The simplest way 

to implement this renormalization of the shift is by sub- 
tracting the equivalent expression for a transparent slab 
with n = 1, 



5E t 



AE* 



A£ trav ( ?i = l) 



(10) 



We decide to place the atom at a position zq > L/2 to the 
right of the slab, substitute the mode functions (|A7IAllj) . 
and get 



d' j k 



En + ui 



I AM 



x eX(k+)eS*(k-)(j2A e 2ik > z ° + 



RX e 



-2ik z ZQ 



(ID 



Similarly, the contribution from the discrete set of 
trapped modes reads 



j^i A, a k z 

(|f k s ACT (ro)| 2 + |f k A ACT (r )| 2 ) |/v| 2 , 



■ u> 

2 



(12) 



which can be written in a more explicit form by substi- 
tuting the trapped modes to the right of the slab from 
Eq. l[Al2l) . 



^ trap = -^E£E 



j^i A, a k z 



d 2 k„ 



x eX(k+)eT(k") \M X \ 



Eji + uj 

2 



L 



S,A 



-2kzo 



(13) 



We note that renormalization makes no difference to the 
trapped-modes contribution to the shift, as the trapped 
modes vanish in the limit n — > 1. 

The total energy shift is obtained by combining the 
travelling-mode contribution Eqs. pip and the trapped- 
mode contribution (fTB")) . At first sight, this is very com- 
plicated, since the former is given by an integral over 
k z while the latter involves a sum over discrete values 
of k z = in that are solutions of the dispersion relations 
(IA19[) . In addition, the shift (fTTj) due to travelling modes 
and its counterpart (|13p due to trapped modes diverge 
when evaluated each on their own, as observed before in 
similar circumstances [l2j]. What helps, is the observa- 
tion that the reflection coefficients R\ have poles in the 
complex k z plane at exactly the values of k z — in that 
are solutions of the dispersion relations (|A19|) . Further- 
more, the residues around those poles are such that the 
sum over k z in Eq. (| 13|) can be re-written as a contour 
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FIG. 2: By closing the contour C one can choose a more suitable 
integration path C In the figure, T = v 'n 2 — 1/n. 



and, in turn. 



Ill = 



C Eji + W 



R T E(k z ,k {l ) e 2ik ^> 
,c kl 



dk z 

dkz T7 I ' 

c Eji + ijj k z 
k 2 

dk z - u -i R TM {k z , fc||) e 2ife ^°. (18) 



i?TM(fc z ,fci|) e 2ifc ^°, (17) 



Before proceeding with the evaluation of Eqs. (jTTJ) and 
f|18[) . we note that both of them are written in terms of 
the position Zq of the atom, measured from the center 
of the slab. In practice, one would of course want to 
know the energy shift of the atom as a function of its 
distance to the surface of the slab, Z = zq — L/2. Writing 
Eqs. (JTTJ) and (fT5| in terms of the atom-surface distance 
Z gives rise to a phase factor e lkzL which we absorb in 
the reflection coefficients (|A8|) by redefining 



integral with the same integrand as in Eq. (fTTj) . Thus the 
sum of Eqs. (jTTJ) and (TlU)) can be combined into a sin- 
gle contour integral in the complex k z plane. In Ref. [7J 
we have shown this to be the case in connection with 
a proof of the completeness of the electromagnetic field 
modes around a dielectric slab, and we refer the reader 
there for details. Using this method to add Eqs. (jTTJ) and 
(fl~3|) . we obtain for the total energy shift 



5E = - 



1 



2(27r)3e 



EE ^ 



dk z — — 

c Eji 



x|/vr e a x (k + )el*(k-)R x (k z M) e 



2ik z ZQ 



(14) 



with the integration path C as shown in Fig. [2] The 
poles of R\ lie on the imaginary axis between and 
iy/n 2 — 1/cy /n, so that C runs above them. To manipulate 
this expression further, we sum over the two polarizations 
and re-arrange the Cartesian components a — x,y, z into 
parallel and perpendicular parts relative to the surface 
of the slab. The double integral in kii can be simplified 
by transforming into polar coordinates and carrying out 
the integration in the azimuthal angle, so that the total 
energy shift reads 



8E 



]=jtt <r=||,_L 



(15) 



with parallel and perpendicular contributions given, re- 
spectively, by 



Sn = 



dfcii fen in and Sj_ = 



dfcy fcy/j. 
(16) 



R\ 



1 



so that Eq. IT8|) turns into 



c 



Eji + oj k 2 



Rtm e 2lk ' z , 



(19) 



(20) 



and similarly for Eq. (I17p . Since the photon frequency is 

given by w = \J k 2 + k 2 , one can identify branch points 

in the integrand at k z = iifey . We choose to place the 
square-root cuts from k z = ifeii to ioo and from k z = —ik 
to —ioo. Apart from this square-root cut, the integrand 
is analytic in the upper half-plane and for Z > van- 
ishes exponentially on the infinite semicircle in the up- 
per half-plane. Therefore Cauchy's theorem allows us to 
deform the original integration path C into a new path 
C that goes round the square-root cut from ifcii to ioo 
(see Fig. [3]). Identifying the correct sheet on each side of 
the cut by demanding that lu > on the real k z axis and 
re-expressing k z = iq, we can work out the integral along 
the path C and obtain 



2K 



dq 



R 



TM 



,-2qZ 



q 2 — k 



2 ^ji II 



(21) 



The calculation of the parallel contribution 7y, Eq. (|17p . 
runs along exactly the same lines. Substituting Eq. (|2"Tj) 
into Eq. (|16p . we see that the two contributions 5y and 
S± to the energy shift (fT5|) are both double integrals, 
over fcy and over q — —\k z . We choose to make a change 
of variables to u = (q 2 — k 2 ) 1 / 2 / Eji and v = k^/Eji, and 
arrive at 



Sn 



du 



dv 

JO 

(u 2 + v 2 )R 



1 



y/u 2 + V 2 1 + ' 



(22) 



TM 



u z R 



TE 



5 



and 



IV. ASYMPTOTIC ANALYSIS 



Si 



dv / du 



1 



Vu 2 + v 2 1 + u< 



Rtm 



x e 



-22£ jiV / « 2 +» 2 



(23) 



These expressions, together with Eq. (jT5J), give us a gen- 
eral formula for the energy shift of a ground-state atom 
in front of a dielectric slab. The shift depends on the 
matrix elements of the atomic dipole between the initial 
state \i) and other states \j) that are coupled to \i) by 
strong dipole transitions. In practice, the sum in Eq. (|15p 
over intermediate states \j) is in most cases dominated 
by a single close-lying state with a strong dipole transi- 
tion to the initial atomic state Alternatively, we can 
use the identity 



I01m ct N)I 2 = ^?I01p.N)I 2 , 



(24) 



where a = e 2 /47reo is the fine structure constant, and 
re- write the energy shift in terms of the matrix elements 
of the momentum operator, 



SE 



2a 



irra 



3 ff=||,± 



EjiSa\v<r\ 



(25) 



In this form the shift is very similar to that of an atom 
in front of a dielectric half-space (cf. Eqs. (2.12), (2.25), 
and (2.26) of Ref. [l(|), except for the different reflection 
coefficients in each situation. 

In order to further analyse or calculate the energy shift 
numerically, it is convenient to transform the double in- 
tegrals in Eqs. ([22]) and ([23]) into polar coordinates by 
substituting u = s cos (f> and v = s sin (j), and then replace 
4> by t = cos 4>. This gives 



a-i 



ds / dt 



, 2 t 2 



— (Rtm 



— t 2 R,TE ) 



1 



/•OO />1 

5x = - / ds / dt 
^ Jo Jo 



1 



(l — t 2 ) Rtm 



with reflection coefficients 

-(n 2 - l)t 2 



Rte = 



R 



TM 



2 + (n 2 - l)t 2 + 2y/l + (n 2 - l)i 2 cothA 
n 4 - 1 - (n 2 - l)i 2 



(26) 



(27) 



(28) 



? 4 + 1 + (n 2 - l)t 2 + 2n 2 s/l + (n 2 - l)t 2 cothA 



and the abbreviation A = LEjiSWl + in 2 — l)t 2 ). Thus 
the energy shift of the atom in front of a dielectric slab 
is given by Eqs. US]) or (J25|, with Eqs. (|26 ]) -([29 ]) . In 
this form the shift is readily computed numerically, as 
we shall do in Section |V] However, to extract important 
physics and be in the position to make quick estimates, 
one should investigate the asymptotic behaviour of the 
shift in various physically significant regimes, which we 
shall do first. 



The nature of the interaction of the atom with the 
slab depends on the separation between them: for large 
separations the interaction is manifestly retarded, but 
for small separations the retardation can be neglected 
and the interaction can be assumed to take place instan- 
taneously. The scale on which one makes this distinc- 
tion of the atom-surface separation being small or large, 
comes from comparing the time 2Z/c that a virtual pho- 
ton takes for a round-trip between atom and surface to 
the time scale of internal evolution of the atom. For the 
atom in state \i) with a strong dipole transition into a 
close-lying state \j), the time scale of the atom's inter- 
nal dynamics is given by %/E^. The ratio of the two 
time scales is 2ZEji in natural units, which can there- 
fore be used as the criterion for retardation: the interac- 
tion is manifestly retarded for 2ZEji ^> 1, as the atom 
has evolved appreciably by the time the virtual photon 
has completed its round-trip, and it can be considered 
non-retarded for 2ZEji <C 1, because the atomic state 
hardly changes while the photon travels to the surface 
and back. In terms of length scales, it is the relative sizes 
of the distance Z of the atom from the surface and the 
wavelength 1 /Eji of the strongest internal transition that 
matter. However, the thickness of the slab L provides a 
third length scale to consider. We shall now consider the 
various asymptotics limits. 



A. Thick slab L > Z 

For a very thick slab, i.e. in the limit L — > oo, we 
can approximate cothA ~ 1 in Eqs. ([28]) and ([29]) . Then 
expressions ([26]) and ([27]) reduce to what they would be 
for a dielectric half-space [l(J. The energy shift for an 
atom in front of a dielectric half-space has been analysed 
in detail previously in both the retarded and the non- 
retarded limits [121 ]. 



B. Thin slab L < Z 

If the atom-surface separation Z is much larger than 
the slab thickness L then the exponentials in Eqs. ([26]) 
and ([27]) effectively cut off the integral at very small val- 
ues of s, so that the argument of the coth in Eqs. I|28p 
and ([211]) stays very small throughout the whole effective 
range of integration. Thus we can approximate the coth 
by its small-argument expansion, cothA ~ 1/A. This 
leads to significant simplifications in Eqs. I[28p and ([2"5]) 
because the square roots in the denominators drop out, 
and we get: 



Rte - 
Rtm 



~{n 2 - l)t 2 



2+{n 2 - l)t 2 + 2/(sLEji) 
n 4 -l-(n 2 - l)t 2 
n 4 + 1 + (n 2 - l)i 2 + 2n 2 /(sLE ji ) 



(30) 
(31) 
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Substituting these into Eqs. (|26|) and (p7|) . we can now 
carry out the t integral. While elementary, this integra- 
tion gives an unwieldy combination of rational functions, 
square roots, and arctan, so that we dispense with writ- 
ing it down. The subsequent integration over s cannot 
be performed analytically, unless we make further ap- 
proximations, which we shall do in the following for the 
retarded and non-retarded limits. 



C. Retarded regime {2ZEji > 1) for a "thin" slab 

If 2ZEji ^> 1 then we can apply Watson's lemma to 
the s integrals in Eqs. (|26]) and (|27|) . So, we substitute the 
approximated reflection coefficients (|50|) and (f3"Tj) . carry 
out the t integration, and then expand the integrand of 
the s integral around s = 0, after which the s integral 
over the leading term becomes elementary. In such a 
way we find that the energy-level shift is given by 



SE 



(n 2 -l)L x 
160ir 2 e n 2 Z 5 ^ 



(32) 

We would like to note that this result is valid for 2ZEji ^> 
1 and Z ^> L, but other than that for any slab thickness 
L. In particular, there is no restriction on LEji, which 
can have any size > or < 1 provided it is much smaller 
than ZEji. In this sense the notion of a "thin" slab is 
slightly misleading in the retarded regime, as any slab of 
finite thickness can be considered thin for large enough 
Z. 

Another interesting aspect of this result is that it shows 
that there is absolutely nothing unusual or non-analytic 
about the limit L — ► 0. A calculation of the Casimir- 
Polder force using field-theoretical means and four pho- 
ton polarizations in a Gupta-Bleuler quantization scheme 
[3] has found different results for different ways of imple- 
menting the boundary conditions on the photon field, and 
the tentative explanation for this discrepancy, as given in 
Ref. [3], has been that these different results should ap- 
ply to thick and thin slabs. However, this explanation 
is inconsistent with our explicit results for the Casimir- 
Polder energy shift for slabs of arbitrary finite thickness. 



D. Non- retarded regime (2ZEji <C 1) 

In the non-retarded limit the interactions between the 
atom and the slab can be approximated as instantaneous, 
and the energy shift in this regime can be calculated by 
considering the limit Eji — > 0. One could take this limit 
in Eqs. flU-flM]), with Eqs. (15]) or (gU), but the calcu- 
lation is much shorter if instead we go back to Eqs. (|16p - 
(TT5)) . which was before we had deformed the contour C 
in the complex k z plane. In the limit Eji — > 0, we get 
uj/(E 31 +uj) — > 1 in Eqs. (TT7|) and (fTgj). so that the square 

root cut due to w = A / k 2 + k 2 disappears from them. 



Instead, we get poles at k z = ±ifc|i . When we close the 
contour C in the upper half-plane, we pick up the residues 
of the integrands at k z = iku, so that Eqs. (fl~7f and (fl8|) 
turn into 



(5 + 9n 2 )|^| 2 + 2(4 + 5n 2 )| M J ; 



2tti Res 



27ri Res 



k 2 RTM{k z ,k\\) 
(k z -ik\\)(k z +ifcy) 

k 2 R T M(k z ,k«) 



(33) 



k z — i/c it 



(k z — ifcii ){k z + ifcii ) 



2ik z Z 



,(34) 



which are straightforward to determine. Substituting the 
results into Eqs. (fil)|) and (TT5)) . we obtain for the energy 
shift in the non-retarded regime 



n 2 - 1 \ - , l9 



167T£() n 2 + 1 

3T1 



x / &kk 2 e~ 2Zk — 
Jo i 



-2kL 



(35) 



The same result could be achieved from purely electro- 
static considerations. The atom can be viewed as a 
dipole, and the dielectric slab can be modelled as a se- 
ries of image dipoles. The energy shift is then just the 
Coulomb interaction energy of the atomic dipole and its 
images. We show in Appendix[B]that an electrostatic cal- 
culation of this sort indeed reproduces the energy shift 

Finally we consider the limit L <^ Z and obtain for 
the non-retarded energy shift of an atom near a thin slab 



3{n 4 



1) L 



2567r£ n 2 Z A 



(36) 



V. SUMMARY AND CONCLUSIONS 

We have obtained a general formula for the energy- 
level shift in a ground-state atom near a non-dispersive 
dielectric slab of refractive index n: Eq. (fT5)) . or alterna- 
tively Eq. (|25p . with the parallel and perpendicular con- 
tributions 5|| and Sj_ given by (f2"6"| and (f2"T|) . respectively. 
While given only as a double integral, it is nevertheless 
in a form that is readily amenable to both numerical cal- 
culations and analytic approximations. We have given 
appropriate asymptotic formulae in both the retarded 
(2ZEji ^> 1) and the non-retarded regimes (2ZEji <C 1). 
For the latter we showed that the result can be repro- 
duced by means of a classical electrostatic treatment. For 
thin slabs the electrostatic energy shift varies as L/Z 4 , 
as shown in Eq. (|3l)j) . 

In the retarded regime, on the other hand, our gen- 
eral formula reduces to Eq. (|3"2"|) , showing that the shift 
behaves as L/Z 5 , provided L ^ Z. For this case, it 
is possible to compare our result with the one given in 
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Eq. (218) of Ref. [g] for the interaction energy between 
a ground-state atom and a magnetodielectric plate, 



U(Z) = 



hca(O) L 

1607T 2 £ Z^ 



14e 2 (0) - 

6(0) 



6^ 2 (0) - 1 
M(0) 



(37) 

where a(0) is the static polarizability of the atom. In or- 
der to compare this result to our result (j32j) . we need to 



substitute e(0) = n for the static dielectric constant and 
/i(0) = 1 for the static magnetic permeability. Further- 
more, the diagonal elements of the atomic polarizability 



E 



2E ji \{jW\i)\' 



v = {x,y,z}, (38) 



so that we get for the static polarizability of the isotropic 
atom considered in Ref. M 



LE.=10 




E 



v = {x,y,z} 



31 



In this language our expression ([32]) reads 



6E = 



a(0) 



1607r 2 n 2 e Z [ 



V - 1)(9 + 14n 2 ), 



(39) 



(40) 



which agrees with Eq. (l37|) upon substitution of e(0) = n 2 
and n(0) = 1. 

The great advantage of our general formulae (TT5]) , (l2l?l) 
and (|27|) is that they make it possible to know how the en- 
ergy shift behaves for various slab thicknesses and values 
of the atom-surface separation Z. Using these formulae 
and standard software packages such as Mathematica or 
Maple, one can easily plot SE for any desired parameter 
ranges. In order to plot some examples in a meaningful 
and informative way, we re-write the energy shift in the 
following form 



slab | 1 2 
H'Z | 



47T£ ^ AttE h Z 4 V II " 11 1 

( 41 ) 

with parallel part and perpendicular contributions de- 
fined by 

Wf ab = MZ A E% 5,| and Wf &h = &AZ A E%S±, (42) 



FIG. 3: The function VF| lab for various thicknesses of the dielectric 
slab, with refractive index n = 2. The uppermost curve is the result 
for a dielectric half-space, W? . 




FIG. 4: The function W| lab for various thicknesses of the slab, 
with refractive index n = 2. The uppermost curve is for a dielectric 
half-space, i.e. for LEji — » oo. 



and the functions Snx given as before in Eqs. (|26p and 
(|2T|) . The motivation for this choice is that (i) W\\ and 
Wx are dimensionless quantities, and (ii) they facilitate 
easy comparison to the standard Casimir-Polder result 
as Wii = 1 = W± for the retarded energy shift of an 
atom in front of a perfect mirror [lH ] . When interpreting 
the plots it is important to bear in mind that one needs 
to multiply with a factor — 1/Z 4 in order to judge the 
distance dependence of the energy shift. For example, the 
functions W^ z (|4"2"1) are linear for small Z, showing that 
the energy shift for small distances behaves as —1/Z 3 , as 
one expects for an electrostatic interaction. 



In Figs. [3] and 0] we have plotted W» 



of ZEji for several slab thicknesses LEji, 



as functions 
while fixing 

the refractive index to n — 2. We have also included 
these functions for the dielectric half-space pUHHi which 
corresponds to the limit LEji — > oo. In Fig. [3] we show 
how the shift varies for different refractive indices if we fix 
the thickness of the slab at LEji = 1. In practice values 
of LEji ~ 10 might be more realistic, but for those the 
energy shift is almost indistinguishable from the one for 
a dielectric half-space, as evident from Figs. [3] and 0J 
Fig. [5] shows W z for various refractive indices n — 
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1.0 




FIG. 5: The function W| lab for a slab of thickness LEji = 1 
and various values of the refractive index, n = 1.5,3,5,10. The 
uppermost curve is the result for a perfect reflector, i.e. for 

n — ♦ co. 



FIG. 7: The function W^ lah as a function of LEji, f° r an atom 
located at a distance ZEji = 8 from the slab and various values for 
the refractive index n = 1.5, 3, 5. In the limit n — * oo this function 
approaches a unit step function. 




o.o -) , 1 , 1 , 1 , 1 , 

2 4 6 8 10 

LE 

I 



FIG. 6: The function Wf ah as a function of LEji, for an atom 
located at various fixed distances ZEji = 1, 5, 10 from the surface. 
The slab has a refractive index of n = 2. 



1.5, 3, 5, 10, and we have also included the limit of a per- 
fect reflector, n — > oo, labeled as W^. Furthermore, one 
can see how the energy shift varies with the thickness of 
the slab LEji. In Fig- El we have plotted W z as a func- 
tion of slab thickness for various fixed surface-atom sep- 
arations, fixing the refractive index at n = 2. For Fig. 
we have fixed the atom's position at ZEji — 8 in the 
retarded regime, and shown how W| lab varies with the 
slab thickness for various values of the refractive index 
n. This shows again that the retarded Casimir-Polder 



force has a well-defined and analytic limit for L — > 0. 

APPENDIX A: MODE FUNCTIONS 

Throughout this paper we have adopted the same no- 
tation as in [7J, where k is the wave vector in vacuum 

k ± = (fe K ,fe J/ ,±fe z ) = (k|| ) ±A; z ), (Al) 

and kj = (ky,±fc z d) is the wave vector inside the di- 
electric slab. The ^-components of the wave vectors in 
free-space and dielectric are related through Snell's law 
by 

fc zd = ^(n2-l)fcf +n2fc2, (A2) 

and in reverse 

k z = i^-^-l)^, (A3) 

which are always positive. 

The vector mode functions can be written as a product 
of a polarization vector and a scalar mode function, 

fkA(r)-e A / kA (r) (A4) 

We work with the transverse electric (TE) mode, for 
which the electric field is perpendicular to the plane of 
incidence, 

e TE = {-A {l )- 1 / 2 (-id y ,id x ,0), (A5) 

and the transverse magnetic (TM) mode, for which the 
magnetic field is perpendicular to the plane of incidence, 

e TM = (AAiir^-a^-cV^Aii) . (A6) 



The momentum space representations (k ) of the po- 
larization vectors are obtained by applying the above dif- 
ferential operators to a plane wave e lk r . 

The scalar mode functions for travelling left-incident 
modes read 



/kA(r) = 



e lk+ ' 
J A e* 
T\e a 



Rxe* 1 
+ Jxe^ 



z < -L/2 

\z\ < L/2 (A7) 

z > L/2 



for any polarization A = TE,TM. The normalization 
constant is N = (2n)~ 3 ' 2 , and the remaining coefficients 
are obtained from the continuity conditions ([3]); in par- 
ticular, 



Rx = 
Tx = 



1 



e 



2ik zl ±L 



rx - 



-ik z L 



^•2 g2i/c Z( -i.Zj 



1-7-2 



r 2 p2ik z aL 
1 A e 



= i(fe z d-fez)i 



(A8) 
(A9) 



where 



fc zd n h z fc zd 

- — and r T M = 2 , , 

TL K z -\- fc zd 



Vzd 



(A10) 



The right-incident modes can be obtained straightfor- 
wardly from the left-incident modes, by simply inverting 
the z-axis and taking z — > — z. 



f£x(*) = N 



T x e^ 
Ixe 1 ^ 



+ Jxe 1 
Rxe ik " 



z < -L/2 
\z\<L/2 (All) 



The trapped modes are given by 



±Lf' A e ik ir r+K2 , 

gikf-r ± e i^ r : 

j S,A iki -r— kz 
X ' 



z > L/2 



z < -L/2 
\z\ < L/2 
z > L/2 



(A12) 



where the ± signs apply to the symmetric (S) and an- 
tisymmetric (A) modes, respectively, and k = |ifc 2 | > 0. 
Note that for trapped modes the polarization vector (|A6|) 
is complex and no longer of unit length. The normaliza- 
tion constants are 



M- 



1 



TE 



Wn 2 f + M? 



Mtm = 



1 



(A13) 
(A14) 



k kf.+n 2 K 2 



and by imposing the continuity conditions ^ we obtain 



the coefficients L 



L, TE 



S.A 



k z dL \ kL / 2 



T A — 



2 cos 



2 e KL / 2 



2n cos 



T A 



= 2ni sin 



K.dL \ kL / 2 



(A15) 
(A16) 
(A17) 
(A18) 



The dispersion relations that arise from the simultane- 
ous application of all matching conditions in Eq. ([3]) to 
the symmetric (S) and antisymmetric (A) modes, with 
two polarizations A each, read 



fc zd tan(& zd L/2) 
-fc zd cot(fc zd L/2) 
-fc zd cot(fc zd L/2)/n 2 
fc zd tan(fc zd L/2)/n 2 



for (S), A = TE, 
for (A), A = TE, 
for (S), A = TM, 
for (A), A = TM, 



where 



- x l(n*-l)kl-kl d 



(A19) 



(A20) 



APPENDIX B: ELECTROSTATIC 
CALCULATION OF THE ELECTROSTATIC 
SHIFT 



In order to have an independent check of our general 
formula for the energy-shift, which in the non-retarded 
limit takes the form (f35|) . we shall derive the same non- 
retarded shift purely by means of a classical electrostat- 
ics. If retardation can be ignored, the energy shift of 
the atom is simply the electrostatic energy of the atomic 
dipole when placed near the dielectric slab. 

If the electrostatic potential $(r) generated by a unit 
point charge at a position r' is known, then the electro- 
static energy of an atomic dipole located at ro is (cf. e.g. 
[l6| for a more detailed discussion): 



i={x,y,z} 

"(El) 

Here the harmonic function (r, r') is the difference be- 
tween the potential <5>(r) generated by the point charge 
at r' and the potential that would be generated by that 
charge in unbounded space, so as to exclude from AE os 
the (infinite) electrostatic self-energies that do not de- 
pend on the relative position of the dipole and the slab. 
As $ir(r,r') is a solution of the Laplace equation and 
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must vanish for z — > ±00, it must be of the form: 



$ ff (r,r') 



dk-, 



dk y e ik * x+ik y y 



Ci(k||,r') e k « z for z< -L/2 

x <( C 2 (k||,r') e fc ii z +C 3 (k||,r') e- fc ii z for \z\ < L/2 
C 4 (k||,r') e~ k ii z for z> L/2. 

The coefficients Ci_4(kii , r') are easily worked out by ap- 
plying the continuity conditions ((3|) to this electrostatic 
problem. Straightforward manipulations then give 



H 



1 £-1 
47T£ £ + 1 J Q 
1 



dk Jo(kp) e 



-k(z+z'-L) 



-2kL 



-l\2 n -2kL 



for z,z' > L/2 , (B2) 



with p — \J{x — x') 2 + (y — y 1 ) 2 . It is instructive to 
rewrite the denominator as a geometric series 



E 



n=0 



and note that [H Eq. 6.611(1 



£+1 



-2kL 



dk Jo{kp) e 



k(z z ) 



\/p 2 + (i - z') 2 ' 



which reveals that $h (r, r') can be understood as being 
due to a series of image charges generated by repeated 
reflections between the two interfaces of the slab fl7j |. 
However, expression (IB2|) is more useful for calculations; 
substituting it into Eq. (|B1[) gives for the electrostatic 
energy shift 



1 £ - 1 



167T£ £ + 1 



(2M 2 + |mii 



< / dk k e 
Jo 



2 -2Zk 



1 - e 



-, (B3) 



z,— 2kL 



which, upon replacing e = n 2 , is in agreement with 
Eq. . 
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